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ABSTRACT 

The  Kolmogorov  Arnold  and  Moser  (KAM)  theorem  states  that  a  lightly  perturbed  Hamiltonian  system  will  have 
solutions  which  lie  on  a  toms.  Earlier  work  by  the  second  author  has  shown  that  most  Earth  satellite  orbits 
perturbed  by  the  geopotential  lie  on  KAM  tori.  The  problem  then  arises  as  to  how  to  convert  the  current  satellite 
tracking  orbits  to  KAM  tori.  A  KAM  torus  is  characterized  by  three  frequencies  and  three  phase  angles.  The 
frequencies  are  essentially  the  rates  of  change  of  the  mean  anomaly,  the  longitude  of  the  ascending  node,  and  the 
argument  of  perigee.  In  this  paper  we  explore  the  determination  of  these  three  rates  from  the  fitting  of  SGP4  two 
line  element  sets  (TLEs),  and  then  constructing  KAM  tori  with  the  specified  frequencies.  The  success  of  this 
process,  and  an  idea  of  the  residual  errors,  can  then  be  obtained  by  comparing  the  SGP4  predictions  with  the  KAM 
torus  predictions.  Second  order  polynomials  are  fit  to  data  from  TLEs  over  18  months  using  a  least  squares 
technique.  The  first  order  coefficients  are  used  as  the  torus  basis  frequencies  while  the  second  order  terms  are  used 
to  account  for  perturbations  to  the  satellite's  orbit  such  as  air  drag.  Four  cases  are  attempted  using  the  Hubble  Space 
Telescope  and  three  rocket  bodies  as  test  subjects.  A  KAM  torus  with  the  desired  basis  frequencies  is  constructed 
and  used  to  predict  satellite  position.  For  the  final  test  case,  this  shows  an  oscillatory  error  with  an  amplitude  of  less 
than  80  km  over  a  period  of  almost  two  years.  The  authors  speculate  that  this  is  caused  by  periodic  lunar  and  solar 
perturbations,  masked  in  the  SGP4  fits  by  frequent  updates. 

1.  MOTIVATION 

Current  methods  of  orbit  propagation  using  a  general  perturbations  technique  such  as  the  SGP4  model  can  only 
accurately  predict  a  satellite’s  position  for  a  few  days  before  an  update  is  necessary.  In  order  to  obtain  the 
information  necessary  to  update  the  model,  a  series  of  observations  made  with  radar  or  electro-optical  sensors  is 
necessary.  If  more  accurate  or  long-term  predictions  are  needed,  for  collision  avoidance  for  example,  numerical 
integration  is  used.  While  advances  in  computational  technology  have  greatly  reduced  the  time  necessary  to 
perform  these  integrations,  they  still  take  hours  to  complete.  In  addition,  to  predict  position  one  week  in  advance  for 
instance,  the  orbit  must  be  integrated  through  all  of  the  intermediate  time  steps  to  produce  the  one  needed  predicted 
location.  A  new  method  with  comparable  accuracy  and  less  computational  and  observational  cost  would  be 
beneficial. 


2.  BACKGROUND 

Kolmogorov,  Arnold,  and  Moser  Theorem,  or  KAM  Theorem,  was  originally  developed  in  the  1950s  as  an  attempt 
to  solve  the  classical  three-body  problem.  In  words,  KAM  theorem  states  that  for  an  integrable  system  in  which  the 
momenta  and  forces  are  invariant  -  that  is,  a  Hamiltonian  system  -  subject  to  small  smooth  perturbations  from 
conservative  forces,  many  solutions  for  the  unperturbed  system  are  also  solutions  to  the  perturbed  system  with  small 
changes.  The  perturbed  system  contains  one  or  more  action-angle  variables  such  that  the  Hamiltonian  can  be 
written  as 


H  (J,<0)  =  hV)  +  efU,<0)  (1) 

Where  /and  co  are  the  action-angle  variables,  A  is  the  unperturbed  Hamiltonian,  /  is  the  perturbing  function  and  e  is 
the  small  perturbing  parameter.  Solutions  to  the  system  can  then  be  shown  to  have  the  characteristic  of  returning  to 
their  initial  position  if  one  angular  coordinate  is  incremented  by  an  integer  multiple  of  some  characteristic  angular 
period  while  the  others  are  held  fixed.  This  type  of  solution  is  defined  as  a  torus,  with  characteristic  frequencies 
defined  by  the  constant  rates  of  change  of  the  angular  coordinates. 
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Previous  work  by  Wiesel  has  shown  that  an  Earth-orbiting  satellites’  orbit  can  be  described  as  a  three-dimensional 
torus  in  six -dimensional  phase  space  [1],  This  is  possible  by  expressing  the  satellite’s  position  in  the  Earth- 
Centered -Earth-Fixed  (ECEF)  reference  frame.  In  this  frame,  the  satellite  position  vector  is  dependent  on  three 
periodic  terms  which  can  be  used  as  approximations  to  the  three  characteristic  frequencies  of  the  torus,  o)t: 


co,  =  M  « 

■y  a 3 


372R© 


2a2(l 


(x)- 


—  fl  —  60  0  ~  — 


3Vm7  2R2 


e2y.s 

2 
© 


(-sin2(i)  -  l) 
cos(i)  —  £U@ 


CO  z,  —  CO 


2a35(l  —  e2)2 

®  (^sin2(0-l) 


(2) 


3  Vm/2^©  c5 

2a35(l  —  e2) 


Where  co©  and  R@  are  the  rotational  rate  and  radius  of  the  Earth  respectively,  /2  is  the  dominant  coefficient  in  the 
geopotential,  //  is  the  gravitational  parameter  and  M,  Q,  co,  a,  e  and  fare  the  familiar  classical  orbital  elements. 

These  frequencies  can  then  be  used  to  form  the  action-angle  variables  for  the  system.  Specifically,  the  angular 
coordinates  are  the  integrals  of  the  three  frequencies  and  are  all  linear  with  time  -  M,  H  —  <u@t,  and  co  -  while  the 
conjugate  momenta  can  be  calculated  using  the  principle  of  least  action  and  have  been  shown  to  be  approximately 
equal  to  the  Delaunay  momenta  [2],  Using  these  action-angle  variables,  the  approximate  Hamiltonian  can  be  written 
as 
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Where  P,  are  the  Delaunay  momenta. 

The  method  of  producing  a  torus  capable  of  accurately  modeling  a  satellite’s  orbit  has  been  developed  and  reported 
on  previously  by  Wiesel  [  1]  [2]  and  can  be  broken  into  3  steps: 

1. )  Numerically  integrate  the  orbit  to  produce  a  dense  set  of  position  vector  data  in  the  ECEF  frame 

2. )  Identify  the  orbit’s  fundamental  frequencies  (the  torus  basis  frequencies)  to  a  high  degree  of  accuracy  using 

a  modified  Laskar  frequency  algorithm 

3. )  Write  each  physical  coordinate  (x,y,z)  as  a  Fourier  series  in  the  three  torus  basis  frequencies 

This  method  has  been  shown  by  Bordner  to  predict  a  satellite’s  position  to  an  accuracy  of  tens  of  meters  over  a 
period  of  six  months  for  certain  types  of  orbits  [3],  Other  research  by  Derbis  and  Little  has  demonstrated  limitations 
of  the  method  in  that  it  cannot  readily  accommodate  non-conservative  forces  such  as  air  drag  or  perturbations  due  to 
station-keeping  maneuvers  on  operational  satellites  [4]  [5], 

The  current  work  builds  on  the  torus  construction  method  described  above  by  using  observational  data  in  the  form  of 
TLEs  to  identify  the  torus  basis  frequencies  and  adding  a  second-order  component  to  the  torus  coordinates  to 
account  for  air-drag. 


3.  METHOD 

For  the  current  effort,  four  test  objects  were  chosen  from  the  satellite  catalog  based  on  the  criteria  that  they  were 
uncontrolled  satellites  in  LEO  orbiting  above  300km  such  that  air  drag  effects  would  be  minimal.  Tab.  1  gives 
general  information  on  the  approximate  orbits  of  each  test  object. 


Tab.  1.  Orbital  information  for  test  satellites 


Name 

Catalog# 

Launch  Date 

Period 

[min] 

Inclination 

[deg] 

Apogee 

[km] 

Perigee 

[km] 

Eccentricity 

Hubble  Space  Telescope 

20580 

4/24/1990 

95.93 

28.47 

566 

561 

3.836E-04 

Thor  Ablestar  Rocket  Body 

59 

10/4/1960 

106.44 

28.25 

1203 

921 

1.893E-02 

Delta  1  Rocket  Body#l 

341 

7/10/1962 

157.52 

44.77 

5619 

949 

2.414E-01 

Delta  1  Rocket  Body  #2 

8133 

8/27/1975 

95.21 

25.30 

700 

357 

2.753E-02 

Approximately  18  months  of  TLE  data  was  collected  for  each  test  subject  from  the  website  http://celestrak.com.  M, 
H,  and  w  data  was  extracted  in  order  to  determine  the  toms  basis  frequencies.  This  was  done  by  plotting  the  data  vs 
time  and  fitting  second-order  polynomials  to  the  data  using  a  least  squares  technique.  The  curve  fits,  and  the 
residuals  from  those  fits  for  the  Thor  Rocket  Body  can  be  seen  in  Fig.  2.  The  first-order  coefficients  from  these 
curve-fits  are  the  desired  torus  basis  frequencies. 

Starting  from  an  initial  position  and  velocity  derived  from  the  initial  TLE  in  the  data  set  and  SGP4  algorithms,  each 
test  subject’s  orbit  was  numerically  integrated  for  a  total  of  one  year’s  time.  The  numerical  integrator  used  takes 
into  account  the  first  twenty  terms  in  the  Earth’s  geopotential  from  the  EGM96  model.  Error  checking  was 
accomplished  by  calculating  the  Hamiltonian  given  in  Eq.  3  at  each  time  step  and  ensuring  its  value  stayed  constant. 
For  all  cases,  error  in  the  Hamiltonian  did  not  exceed  10"13. 

The  numerical  data  was  then  analyzed  to  identify  the  orbit’s  fundamental  frequencies  using  the  modified  Laskar 
frequency  algorithm  developed  and  applied  previously  by  Wiesel  and  Bordner  [1][3],  Fig.  1  shows  the  partial 
frequency  spectrum  of  the  Thor  Rocket  Body  orbit  and  identifies  the  major  peaks  near  integer  multiples  of  the  Mean 
anomaly  frequency  (coj  and  2oo i). 


Partial  Frequency  Power  Spectrum  for  Thor  Rocket  Body 


Fig.  1.  Partial  frequency  spectrum  for  the  Thor  Rocket  Body 

Having  identified  the  fundamental  frequencies  of  the  numerically  integrated  orbit,  the  numerical  data  was  fit  to  a 
Fourier  series  in  the  three  basis  frequencies.  This  Fourier  series  expresses  physical  coordinates  in  terms  of  torus 
coordinates  and  is  of  the  form 


q  =  Cjcosj  •  Q  +  Sjsinj  •  Q ) 

J 

Where  if  is  a  vector  containing  physical,  Cartesian  coordinates,  /  is  a  vector  summation  index  with  each  vector 
corresponding  to  a  linear  combination  of  the  three  basis  frequencies  identifying  a  peak  in  the  power  spectral  density 
function  plot  of  the  numerically  integrated  trajectory,  C and  S are  matrices  of  Fourier  coefficients  and  Q  is  a  vector 
containing  torus  coordinates. 

The  initial  torus  basis  frequencies  were  then  compared  to  the  desired  torus  basis  frequencies  obtained  from  the 
TLEs.  Following  the  initial  torus  construction,  the  differences  between  the  two  sets  of  frequencies  were  not 
negligible.  In  order  to  correct  this  difference,  a  frequency  matching  iteration  was  performed.  This  iteration  was 
accomplished  by  calculating  the  change  in  initial  conditions  (position  and  velocity)  that  would  give  the  desired 
change  in  basis  frequencies.  For  this  work,  it  was  assumed  that  the  initial  position  derived  from  the  initial  TLE  in 


the  data  set  and  SGP4  algorithms  was  accurate,  and  therefore  held  fixed,  while  the  initial  velocity  was  allowed  to 
change.  The  reference  (first-iteration)  toms  basis  frequencies  are  grouped  into  a  vector,  w,  while  the  desired  torus 
basis  frequencies  from  the  TLE  data  are  grouped  into  a>0.  The  frequency  error  can  then  be  written  as 


dal  _> 

Aa>  —  a)  —  u>n  =  — ^-A  P 
dP 


Where  P  is  the  vector  of  the  reference  toms  momenta  and  A P  is  the  unknown  momentum  offset  between  the 
reference  and  desired  toms.  Eq.  5  can  be  solved  for  A P  which  yields 


.»  dal  1 
AP  =  —  Aal 
dP 
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The  Jacobian,  dal/ dP  can  be  found  analytically  in  terms  of  the  classical  orbital  elements  by  taking  the  torus 
momenta  to  be  the  DeLaunay  momenta  and  writing  the  torus  frequencies  as  partial  derivatives  of  the  system 
Hamiltonian,  K,  given  previously  in  Eq.  3: 


_  dK 
(0  = 

dP 

Therefore,  the  Jacobian  can  be  written 

dal  _  d2K 
dP  ~  dP2 

In  this  way,  a  known  offset  in  torus  basis  frequencies  can  be  expressed  as  an  offsets  in  the  torus  momenta. 

In  order  to  form  the  new  toms  with  the  correct  basis  frequencies  (those  derived  from  the  TLEs),  a  new  numerical 
integration  must  be  carried  out  starting  from  some  new  initial  position  and  velocity  conditions.  Therefore,  the  offset 
in  toms  momenta  must  be  expressed  as  a  change  in  the  physical  position  and/or  velocity.  For  this  effort,  as  noted 
earlier,  the  initial  position  was  held  fixed  while  the  initial  velocity  was  allowed  to  change. 


(7) 


(8) 


Let  the  physical  and  torus  state  vectors  be 


X  —  [rT  vT]T  =  [xy  z  vxvy  vz]T 

y  =  [QT  pt]t  =  [M  n  a>  p1  p2  p3]t 


(9) 


Then,  the  change  in  the  physical  state  vector  can  be  written  in  terms  of  a  change  in  the  torus  state  vector  as 


Where  AT  =  [AQr 


A PT  ■  The  Jacobian  dX/dY  can  be  found  analytically  as 


dX  _  IdY dZ\ 
dY  \dZdXj 
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where 


Z=[Mna>aei]T  (12) 

Is  a  vector  containing  the  COEs.  For  a  detailed  derivation  of  the  content  of  dX/dY  see  [6],  Then,  setting  the 
change  in  initial  position,  A r  —  0,  Eq.  10  can  be  re-written  as  two  linear  equations 


(13) 


Af  =  0  =  AAQ  +  BAP 


Av  —  CAQ  +  DAP 


Where  A,  B,  C  and  D  are  quadrants  of  dX/d  Y: 
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Given  A P  from  Eq.  6,  the  change  in  initial  velocity  needed  to  give  a  desired  change  in  torus  basis  frequencies  can  be 
found  by  first  solving  Eq.  13  for  A Q  and  substituting  into  Eq.  14  which  gives 


Av  —  (— CA-1#  +  D)AP  (16) 

Having  now  obtained  a  new  initial  physical  state,  X0,  the  torus  construction  process,  starting  with  a  new  numerical 
integration,  can  be  repeated  to  find  a  new  torus  with  updated  basis  frequencies.  This  algorithm  was  repeated  until 
the  maximum  basis  frequency  error  was  down  to  10"12  rad/TU.  With  an  updated  torus  with  basis  frequencies 
matching  those  obtained  from  the  TLE  curve  fits,  the  satellite’s  physical  position  was  extracted  using  the  torus 
Fourier  series  of  Eq.  4  and  expressing  the  torus  coordinates  as  the  quadratic  functions  of  time  found  in  the  TLE 
curve  fits. 


4.  RESULTS 

The  process  described  above  was  attempted  for  each  of  the  four  test  cases  from  Tab.  1.  For  both  Delta  Rocket 
Bodies,  curves  were  fit  to  the  TLE  data  set,  however  the  residuals  between  the  data  and  the  curve  fits  for  Mean 
anomaly  were  as  high  was  2  radians  for  Delta  Rocket  Body  #2  and  as  high  as  30  radians  for  Delta  Rocket  Body  #1. 
For  Delta  Rocket  Body  #1,  the  authors  believe  the  poor  fit  is  a  result  of  the  way  TLE  data  is  formatted.  Specifically, 
both  revolution  number  and  epoch  time  are  based  off  of  ascending  node  crossing.  This  can  lead  to  errors  in  the 
revolution  number  count  if  the  satellite  is  off  slightly  from  the  ascending  node.  For  Delta  Rocket  Body  #2,  the  poor 
fit  is  believed  to  be  due  to  the  relatively  low  orbit  (perigee  at  357  km),  the  approximately  350  km  difference  between 
perigee  and  apogee,  and  the  high  area-to-mass  ratio  of  an  empty  rocket  body  resulting  in  significant  air  drag  effects 
that  vary  widely  throughout  the  orbit. 

The  Hubble  Space  Telescope  (HST)  is  in  a  nearly  circular  orbit  with  an  eccentricity  of  approximately  0.0004  at  an 
altitude  of  approximately  560  km.  It  was  expected  then  that  the  air  drag  effects  experienced  by  the  HST  would  be 
less  pronounced  and  nearly  uniform  throughout  the  orbit.  This  was  the  case,  as  this  time,  the  curve-fits  left 
maximum  residuals  of  only  approximately  0.1  rad.  With  good  curve-fits,  and  therefore  reliable  basis  frequencies, 
the  torus  construction  process  was  attempted,  however  the  algorithm  to  identify  the  basis  frequencies  of  the 
numerically  integrated  trajectory  failed  to  correctly  identify  the  smallest  basis  frequency,  co3.  The  authors  believe  the 
cause  of  this  difficulty  to  be  the  low  eccentricity  of  the  HST’s  orbit.  In  the  torus  construction  algorithm,  small 
eccentricities  lead  to  singularities,  similar  to  those  encountered  with  classical  perturbation  theory.  This  problem  has 
been  encountered  in  earlier  research  by  Bordner  [3]. 

The  Thor  Rocket  Body  test  case  had  a  much  larger  eccentricity  than  the  HST  (0.02)  and  was  at  a  higher  altitude  than 
the  Delta  Rocket  Bodies.  As  a  result,  the  problems  in  curve-fitting  and  torus  construction  seen  in  the  previous  cases 
were  not  encountered.  Fig.  2  shows  the  curve  fits  for  the  Thor  Rocket  Body;  this  time,  maximum  residuals  were  on 
the  order  of  10"3  rad.  From  these  curve  fits,  desired  torus  Basis  frequencies  were  determined  and  are  shown  in 
Tab.  2. 


Tab.  2.  Thor  Rocket  Body  desired  torus  basis  frequencies 


col  [rad/TU] 

co2  [rad/TU] 

<o3  [rad/TU] 

7. 9  37272207653  70  E-Ol 

-5. 9670592 1687582  E-02 

1.3  675032 1345 160E-03 

A  total  of  four  iterations  of  the  frequency  matching  algorithm  were  necessary  to  match  the  torus  basis  frequencies  to 
within  the  deshed  tolerance.  Tab.  3  shows  the  progression  of  torus  frequencies  along  with  the  residuals,  calculated 
from  the  desired  TLE  curve-fit  frequencies,  and  Tab.  4  shows  the  initial  velocity  changes  made  in  each  iteration  of 
the  frequency  matching  process.  The  magnitude  of  the  total  velocity  change  was  approximately  0.25  percent  of  the 
initial  velocity  magnitude  obtained  from  the  TLEs  and  SGP4. 


Tab.  3.  Torus  frequency  matching  for  Thor  Rocket  Body 


col  [rad/TU] 

u>2  [rad/TU] 

u)3  [rad/TU] 

Initial  Torus 

7.9373203844 1640E-01 

-5. 9  6704  9  0942  2 124  E-02 

1.3  6749993  5865  27  E-03 

Iteration  1  Torus 

7.937 1878014631 2  E-01 

-5.96705  786750 190  E-02 

1.36747666904746  E-03 

Iteration  2  Torus 

7.93727 1923  6 1304E-01 

-5.96705  9  24650 170  E-02 

1.367503  86446849  E-03 

Iteration  3  Torus 

7.93727  2  2 1037809  E-01 

-5.9670592 1734135  E-02 

1.367503  2  2508483  E-03 

Iteration  4  Torus 

7.937272  2077 1728  E-01 

-5.96705  9  2 16877  51  E-02 

1.367503  2 1348700  E-03 

uil  residual  [rad/TU] 

u)2  residual  [rad/TU] 

0)3  residual  [rad/TU] 

Iteration  1  Torus 

8.4406190579 1449  E-06 

-1.34937392023  970  E-08 

2. 6  5  44404 1400 165  E-08 

Iteration  2  Torus 

2.8404065988 141 3  E-08 

2.96258795273729  E-10 

-6. 5 10 168899 19 101  E-10 

Iteration  3  Torus 

-2. 72439071302699  E-10 

4. 655  29698 123  746  E-12 

-1.1 633  2  299696098  E-ll 

Iteration  4  Torus 

-6. 35802521742335  E-12 

1.689620665  60 141 E-14 

-3. 53  9985 145  295  5  2  E-14 

Tab.  4.  Changes  to  initial  velocit; 

y  in  Thor  Rocket  Body  torus  fitting  process 

Vx  [DU/TU] 

Vy  [DU/TU] 

Vz  [DU/TU] 

Initial  Velocity 

-9.980 15949  328921  E-02 

-8. 11 259 14 1324669  E-01 

4. 393293 18458694  E-01 

Iteration  1  AV 

2. 45  9403  8  6  5407  20  E-03 

-2. 5855  642003  2001  E-04 

8. 5 14642  3  948  6448  E-05 

Iteration  2  AV 

-1.0948665  2807646  E-04 

1.50534 1487 14208  E-05 

-3. 3  5  9447  5  846 148  5  E-06 

Iteration  3  AV 

3.859300490769  70  E-06 

5. 547342066593  74  E-07 

1.45063 14 13075  29  E-07 

Iteration  4  AV 

3. 2043  3729724457  E-08 

5. 37472590607683  E-09 

3.03  7792 1205 1191  E-09 

Total  AV 

2. 34602586873582  E-03 

2.42942896228060  E-04 

-8. 1935077  2974792  E-05 

The  final  iteration  with  matched  basis  frequencies  was  then  used  to  predict  the  satellite  location,  and  this  predicted 
location  was  compared  with  the  location  obtained  using  SGP4  algorithms  at  the  epoch  time  of  each  TLE  in  the  data 
set.  The  first  comparison  was  done  calculating  the  torus  coordinates,  Q,-,  as  linear  functions  of  time;  that  is: 


Qi 


0)it  +  &)j0 


The  results  of  this  comparison  can  be  seen  in  Fig.  3. 
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Fig.  2.  TLE  curve-fits  for  the  Thor  Rocket  Body 


Residuals:  SGP4  Position  -  Torus  Position  (1st  Order) 


Fig.  3.  Toms  position  error,  calculated  using  first-order  torus  coordinates 

It  is  evident  that  the  torus  in  this  configuration  provides  a  poor  position  prediction  with  the  position  in  error  growing 
to  a  magnitude  of  nearly  1500km  after  18  months.  This,  however,  is  expected  as  the  TLE  analysis  showed  that  the 
torus  coordinates  increase  as  quadratic  functions  of  time.  To  compensate  for  this,  the  torus  coordinates  were  re¬ 
calculated,  this  time  using  the  entire  quadratic  curve-fit  at  each  time  step.  The  results  of  this  can  be  seen  below  in 
Figs.  4  and  5. 


SGP4  Position  -  Torus  Position  (2nd  Order) 


Fig.  4.  Toms  position  error,  calculated  using  second-order  torus  coordinates 

These  figures  show  that  the  maximum  position  error  magnitude  is  now  only  approximately  78  km,  and  there  is  no 
apparent  linear  or  quadratic  growth  in  the  error  as  was  seen  in  the  previous  comparison  (Fig.  3). 


Magnitude  of  Position  Difference:  SGP4  Position  -  Torus  Position  (2nd  Order) 


Fig.  5.  Magnitude  of  toms  position  error,  calculated  using  second -order  torus  coordinates 

In  order  to  further  understand  the  causes  of  the  remaining  error  in  position  prediction,  a  frequency  analysis  was  done 
on  the  position  error  data.  Figs.  6  and  7  show  the  results  of  this  analysis  and  illustrate  that  the  position  error  in  all 
three  physical  coordinates  oscillates  predominately  at  the  third  torus  basis  frequency  co3,  and  the  X  and  Y 
coordinates  have  a  secondary  oscillation  frequency  equal  to  the  second  torus  basis  frequency,  CO2. 


Fig.  6.  Frequency  spectrum  of  torus  position  error  -  co3  detail 


Fig.  7.  Frequency  spectrum  of  torus  position  error  -  co2  detail 

5.  CONCLUSIONS 

The  results  discussed  above  show  a  number  of  promising  results.  First,  the  HST  and  Thor  Rocket  Body  test  cases 
show  that  it  is  possible  to  extract  accurate  KAM  torus  basis  frequencies  from  observational  data  in  the  form  of 
TLEs.  Second,  a  method  of  making  slight  changes  to  torus  basis  frequencies  was  demonstrated  by  equating  a 
desired  change  in  basis  frequencies  to  a  small  change  in  initial  velocity.  This  capability  allows  the  current  torus 
construction  method  to  be  used  in  tandem  with  the  TLE  frequency  identification  routine  to  obtain  a  torus  with  the 
desired  basis  frequencies  for  a  given  satellite.  The  difficulties  encountered  with  TLE  curve  fitting  for  the  two  Delta 
Rocket  Bodies  show  that  there  may  be  limitations  in  the  orbital  regimes  for  which  the  current  method  of  basis 
frequency  extraction  from  TLEs  can  be  useful;  specifically,  orbits  in  which  air  drag  effects  are  moderate  to  minimal 
with  little  variation  throughout  the  orbit. 

It  is  the  opinion  of  the  authors  that  the  residual  position  error  seen  in  the  Thor  Rocket  Body  test  case  is  due  to  two 
sources.  First,  the  frequency  analysis  of  the  error  discussed  above  shows  that  the  torus  construction  method  may  not 
be  adequately  characterizing  the  contribution  of  the  two  smallest  basis  frequencies  to  the  satellite’s  position. 

Second,  perturbations  from  the  Sun  and  Moon  are  currently  not  included  in  the  torus  model.  These  perturbations  are 
masked  in  SGP4  and  the  TLEs  due  to  frequent  updates.  Current  research  is  being  done  by  Wiesel  to  address  both  of 
these  sources  of  error  by  further  refining  the  torus  construction  algorithms  to  more  accurately  model  the  effects  of 
co2  and  co3  and  by  incorporating  Solar  and  Lunar  effects  into  the  model. 
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